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ESTIMATE OF THE PHYSICAL 
EXTENT OF THERMAL PLUMES 
AT LAKEVIEW 

ABSTRACT 

The cooling water plumes from an operating thermal 
electric generating station at Lakeview, west of Toronto, 
on Lake Ontario, were defined by water temperature measure- 
ments. The generating station discharged between 1400 to 
2500 cfs of cooling water at 15°F above intake water tempera- 
ture. The temperature measurements were made throughout the 
year at various depths in proximity of the discharge to define 
the physical cooling water plume in the lake. A method for 
estimating temperature isotherms was developed by statistically 
analyzing the temperature measurements. Equations were developed 
which defined plume area as a function of distance from the 
discharge, temperature as a function of plume area and plume 
centreline as a function of wind. Different relationships were 
developed for summer and winter conditions based upon the ambient 
lake temperature. The most significant cooling resulted from 
mixing with lake water. 



ESTIMATE OF THE PHYSICAL 
EXTENT OF THERMAL PLUMES 
AT LAKEVIEW 

INTRODUCTION 

With the advent of large thermal electrical 
generating stations employing Great Lakes water for cooling, 
there has been increasing concern on the possible detrimental 
effects from the addition of heat to the water environment. 
Research on the effect of heat inputs on the water environ- 
ment has been extensive (Parker , 1969 and ASCE, 1967), 
However, the methods for predicting the physical extent of 
cooling water plumes on lakes are notably few. This aspect 
is considered essential by the Ministry of the Environment. 
Once defined, the local effects of the heated plume on the 
water environment could be examined. 

Studies were undertaken by the Hydro-Electric Power 
Commission of Ontario (H.E.P.C.) to define the physical extent 
of the cooling water plume at the Lakeview generating station 
on Lake Ontario. Lakeview discharged approximately 1,400-2,500 
cfs at IBT above intake water temperature. Initially, the 
studies consisted of measuring water temperatures with a 
telethermometer at predetermined locations in proximity of 
the discharge throughout the year. Normally, a survey took 
less than 3 hours. Meteorological data collected during the 



surveys consisted of wind and air temperature determinations. 

This report discusses the development of a first 
order estimate of the spatial distribution of isotherms 
within the cooling water pliome. The estimate is a first 
order approximation based on limited survey data. A better 
estimate could be made by solving the heat budget equation. 
However, the solution of the heat budget equation requires 
more extensive measurements of wind, relative humidity, air 
and water temperatures, net solar radiation and water movement 
patterns. Measurements of these parameters are now being made 

by H.E.P.C. 

Csanady (1970) recently showed that rate of heat 
lost to atmosphere is one-third of that lost to the lake in 
studies at Douglas Point on Lake Huron. Palmer (1969) also 
showed that excess heat in a plume is dissipated through hori- 
zontal mixing with the main lake. It will be shown that the 
principal cooling mechanism appears to be mixing with the 
ambient lake water. These results are in conflict with Brady 
(1969) and Parker (1969) who both indicate that the main cooling 
is atmospheric. Eiesbol (1971) in studies for a cooling water 
discharge on Dardanelle Lake in Arkansas found that cooling by 
entrainment of lake water is the main mechanism in the first 
several thousand feet while atmospheric cooling prevails farther 
away from the discharge point. 



EXPERIMENTAL METHOD 

Water temperature surveys were conducted by H.E.P.C. 
with a small fast boat on relatively calm days. Surface 
and some depth temperatures were measured at the intersection 
points of a 1,000 feet square geometric grid. The intersection 
points were defined by visual alignment on shore range poles 
with some assistance from shore based optical transits when 
required. All temperatures were measured with a Yellow Springs 
thermister telethermometer having an accuracy of 0.2"*F. 

An average water temperature survey was completed 
in less than three hours. During the period of survey, the 
following parameters other than water temperature were measured 
several times on shore and then averaged. 

1. Air temperature in °F 

2. Wind speed in mph 

3. Wind direction in degrees from north 

4. Rate of discharge from outfall in cfs. 

Observed isotherms were plotted for each survey. The centre- 
line of the plume was drawn by joining the farthest points of 
each isotherm from the outfall. The area enclosed by each 
isotherm was measured by planimeter. 

PRELIMINARY INVESTIGATIONS 

Several different methods for predicting lake tempera- 
ture isotherms from the available data were tried with varying 



degrees of success in terms of matching actual observed iso- 
therms. Some of the methods tried are described briefly. 

Gifford's Normal Distribution Dispersion Model 

This method assumes that the concentration distri- 
bution of any continuously discharged material in a moving 
stream will be normally distributed across any plane perpen- 
dicular to the moving stream. As the form of the concentration 
distribution is defined, i.e. a normal distribution, it is only 
necessary to define spread (variance) and current as functions 
of downstream location to determine concentrations at different 
distances from the source. This method has been applied success- 
fully in meteorology (Giffordy 2955), ocean discharges (Foxworthy , 
1966) and lakeshore areas (Palmer, 2969). 
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This method was applied with some success to a cooling water 
plume in May 1968 (Falmer, 1969). The method considers 
dilution with ambient lake water. Since lake cooling has 
been shown as the principal cooling mechanism at least in the 
first few thousand feet, this method could be used effectively 
provided current velocities are known. Current velocities in 
near shore areas with shore topographical effects, could not 
be satisfactorily predicted from wind velocities. Therefore, 
the method was rejected. 

John Hopkin's Method 

This method is based on statistical fitting of the 
principle mechanisms responsible for cooling in small ponds 
or lakes (Brady, 2969). The technique was successfully applied 
to experimentally determined plumes with the computed parameters 
comparing favourably with other studies. However, it was neces- 
sary to define the portion of the lake affected by the plume to 
apply the method. For small lakes, this is simply the total 
surface area. As it was not possible to satisfactorily define 
the area of the lake affected, this method was abandoned. 

PROPOSED METHOD 

The above methods described briefly could not be used 
successfully with the available survey data. Therefore, a 
statistical approach was used to develop a prediction model for 



the thermal plume. The data used is listed in Appendix III. 
Equilibrium temperatures (E) were computed for each 
survey utilizing the net radiation absorbed, heat exchange 
coefficient and dew point temperature (Brady, 1969). In the 
absence of measurements of the above parameters during field 
surveys, net radiation absorbed and dew point temperature were 
obtained from the Atmospheric Environment Services and heat 
exchange coefficient interpolated from Brady's (1969) design 
chart. Equilibrium temperatures were then utilized instead 
of ambient lake water temperature to develop statistical 
equations. However, these equations were statistically very 
weak with small multiple correlation coefficients and large 
standard errors. The weaker equations may be due to errors 
in computed equilibrium temperatures. With the absence of 
field measurements for net radiation absorbed, dew point and 
heat exchange coefficient, equilibrium temperature, the 
computation is not considered reliable. Therefore, equations 
with ambient lake temperature were then developed. 

The field survey data were divided into two seasons, 
i.e. winter and summer since the entire data set as well as 
data division into three seasons (winter, summer and spring 
fall) produced statistically weaker equations with large 
standard errors of estimates. The equations were improved 
by separating the data into 'winter' and 'summer' sets. Since 
much of cooling is due to mixing with lake water and air 



temperatures are very variable, ambient lake water tempera- 
ture was chosen as tlae criterion for separation of data. 
Lake temperatures of 42 '^F, 4 5°F and 4 7°F were tried for 
separation of data. 4 ST was found to be the best since it 
produced statistically best prediction ecpiations for both sets. 
Best fit equations were developed for the following: 

a. area (A) contained by isotherms and the 
distance (X) (the farthest point on iso- 
therm from the outfall, measured along 
the plume centreline) 

b. area (A) and temperature difference 
(AT_T. ) between isotherm corresponding 
to area and ambient lake temperature 

e. distance (X) and temperature difference 

From the first two relations, it is possible to 
define plume area and temperature at any distance (X) from the 
outfall. These equations are necessary to determine the plume 
temperature and area at any distance from outfall to enable 
plotting of the plume. The relation between X and AT has 
been obtained for completeness, but cannot be used to determine 
isotherm shapes since it does not define isotherm area. Various 
shapes, namely circular, elliptical and one leaf of multi-leaved 
rose were investigated as possible shapes for the predicted isotherm 
One leaf of multi-leaved rose (Figure 1} was found to approach the 
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observed isotherm shape best and was adopted in the present 
model. The direction of the plume centreline was determined 
from the wind vector. This method will be described in the 
section 'Definition of Plume Centreline'. 

Plume Area and Distance 

Relations have been found for area or spread or 

variance as a function of distance from a source by many 

experimenters (ALsafar, 1966; Csanady^ 1964; Okubo, 1968; 

Foxworthy , 1966 and Murthy , 1969). In the case of nearshore 

areas of the Great Lakes, the relations appear to be general 

in slope for a particular location, but vary in intercept for 

the maximum, mean and minimum conditions (Palmer and Izatt, 

1970). On this basis, it is probable that a statistically 

significant equation could be determined relating area 

(contained by the isotherms) of the plume and distance from 

outfall along the plume centreline. The best fitted equation 

(see Table 1 and Figure 2) for the suiraner data set was: 

Area = -0.696 + 1.384 (distance) (2) 

with F, Qj. = 117.1: multiple correlation 

coefficient = 0.74: standard error = 3.11: 

standard deviation area - 4.62. 

The best winter data equation (see Figure 3) was: 

Area = -4.047 + 2.43 (distance) ... (3) 

with F, ^T = 28 3.6: multiple correlation 
1,63 ^ 

coefficient = 0.90: standard error = 2.49: 
standard deviation area = 5.80. 



It is iiranediately apparent that there was a statistically 
significant relationship between the areas contained by the 
isotherms and the distance from the outfall. The winter 
equation was statistically stronger, when lake temperature 
variations were smallest. 

Isotherm Area and Temperature Difference 
(Plume and Ambient Lake Tanperature) for 
Flows between 1400 and 2500 cfs 

The cooling of the plxime to either the atmosphere 
or the lake must be related to the plume area and rate of 
discharge. The entrainment of lake water at ambient temperature 
is, partly a function of the pliane area in contact with the lake 
water. Similarly, the amount of atmospheric cooling is, in part, 
a function of the surface area. If one assumes that the plume 
area in contact with the atmosphere is nearly the same as that 
in contact with the lake, one would expect a relationship to 
exist between the area contained by an isotherm, rate of dis- 
charge and the temperature difference between plume and ambient 
lake temperature. This seems to be a valid assumption as most 
of the thermal plume appears to be confined to the surface regions 
of the lake. Statistically significant relationships were found 
(see Table 3 and Figure 4). The best equation for the summer 
data including a flow term was: 

Area = 90.53 - 8.86 In (flow x temperature 
difference) ... (4) 
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with F^ Qc - 135.2: multiple correlation 
coefficient = 0.77: standard error = 2.99: 
standard deviation area = 4.62. 
the best winter equation including a flow term (see Figure 5) 
was: 

Area = 74.0 - 6.88 In (flow x temperature 
difference ... (5) 

with F, fro - 126.3: multiple correlation 
coefficient = 0.82: standard error = 3.37: 
standard deviation area = 5.80. 
The above equations may be used if the cooling water discharge 
lies between 1,400 and 2,500 cf s . If the rate of discharge is 
outside this range , equations developed in the next section 
should be used. 

The question immediately arises whether a better 
equation could be determined if temperature differences between 
the plume and atmosphere as well as the temperature difference 
between ambient lake and plume were included. Statistical 
fitting of the summer and winter data including both temperature 
differences is tabulated in Table 2. Significant equations were 
found relating the temperature differences and isotherm area 
and for the summer data the standard error was slightly less. 
However "t" tests on the significance of the two temperature 
differences indicate that the air-plume difference was not 
significant for the summer data and not significant at a = 0.01 
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for the winter data. The predominate cooling mechanism in 
these surveys was lake water cooling. As the method proposed 
here is a first order estimate, it was decided for simplicity 
to ignore atmospheric cooling and utilize the relationship 
between isotherm and lake-plume temperature difference. 

Isotherm Area, Flow and Temperature Difference 
(Plume and Ambient Lake Temperature) for any Flow 

As equations of Table 3 apply when cooling water 
discharge (flow) is 1,400-2,500 cf s , other equations must be 
developed for all rates of flow. Such statistically signifi- 
cant equations {see Table 5) for summer and winter data, were 
then developed, relating Q/A (flow/plume area) to the ATp^^ 
(temperature difference) . The best equation for the summer 
data was : 

In (flow/area) = 1.91 + 2.17 In (temperature 
difference) ••• (6) 

with F, Q^ = 82.0: multiple correlation 
coefficient = 0.68: standard error = 0.96: 
standard deviation (flow/area) = 3.70. 
The best equation for winter data was: 

In (flow/area) = 3.11 + 1.4 2 In (temperature 
difference) ••• (?) 

with F, ,- = 89.0: multiple correlation 
1 , b J 

coefficient = 0.76: standard error = 0.77: 
standard deviation (flow/ area) = 3.30. 
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The above equations (Table 5) should be used only when the 
flow is outside the range of 1,400 - 2,500 cfs. If the flow 
is between 1,400 - 2,500 cfs, the equations of previous 
section (Table 3) are preferred. 

As the dependent and independent variables in 
Tables 3 and 5 are different, comparison of the standard 
errors (in Tables 3 and 5) associated with different vari- 
ables is misleading. 

Isotherm Distance and Temperature Difference 
(Plume and Ambient Lake) 

Although this relationship was not required to 
predict the isotherms, it is useful to obtain rapid estimates 
of maximum distances isotherms will occur away from the dis- 
charge point, without plotting actual isotherms. The equations 
are summarized in Table 4. It is observed that the equations 
were not as statistically significant as the other fits and the 
standard errors over standard deviation of distance were larger 
than the other fits. The best summer equation including a flow 

term is : 

Distance = 33.08 - 3.0 In (flow x temperature 

difference) ••• (^^ 

with F, =29.0: multiple correlation 
coefficient = 0.48: standard error = 2.18: 
standard deviation of distance = 2.48. 
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The best winter equation including a flow term was: 
Distance = 26.88 - 2.30 In (flow x 
temperature difference) ... (9) 

with F, c^ = 73.4: multiple correlation 
coefficient = 0.73: standard error =1.48: 
standard deviation = 2.16. 

The winter equation is the significantly better prediction 

equation (Figure 6). 

Definition of Plume Centreline 

The curvature or radius of the plume was defined by 
fitting the measured radius of the plume centreline and the 
absolute value of the resolved component of the wind acting 
parallel to the shoreline. The locii for the centres for 
the plume radii were determined by trial-and-error and appear 
in Figure 7. The generality of these locii in terms of other 
locations is not known particularly as there was an effect on 
the plume from the cooling water intake channel configuration. 
The best fit relationship for radii and wind component appears 
in Figure 8, 

R = 500 jWp| ... (10) 

where R = plume centreline radius in ft. 

Wp = wind component along the shore in mph. 
It is apparent that the point scatter is large for equation (10) 
The standard error of estimate can be improved by defining a 
unique locus of the centre of the radius for each experiment. 
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However, this is of little value for a technique of represen- 
ting the thermal plume. Other methods involving current 
measursnents could not be used as the measurements are in- 
complete. 

Step Description of Method 

From the meteorological data, for the period under 
consideration, the wind component along the shoreline was 
computed and thence the radius of plume centreline from the 
equation (10). The plume centreline is drawn according to 
the definition of Figure 7. 

For the existing lake temperature, falling in 'summer" 
or 'winter' regimes and the flow rate, the appropriate equations 
were selected from Tables 1, 3 or 5. Plume temperatures and 
areas were computed at various distances from outfall and along 
the plume centreline for a particular discharge. For any distance, 
the area computed represented the plume contained by the corres- 
ponding temperature isotherm. One leaf of the multiple- leaved 
rose was then fitted to the distance and plume area. A step by 
step computation appears in Appendix I and II. 

Three examples of the method are presented in Figures 
9, 10 and 11 for lake temperatures of 33°F, 46.5°F and 71.7°F 
and wind vectors acting in three different directions. The 
example of Figure 11 is repeated, using equations of Table 5 
instead of Table 3 and presented in Figure 14 (see Appendix II) . 
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The actual field survey results appear as solid line iso- 
therms whereas the predicted plumes are broken lines. The 
shaded area for the plume predictions represent one standard 
error of estimate for the plume area with isotherm enclosing 
the smaller area representing the mean value. Similarly, 
the temperature range is one standard error of estimate with 
the upper value representing the mean value. Errors of 
locating the plume centreline have not been included and the 
centreline represents the mean value location. 

DISCUSSION 

The isotherm predictions in Figures 9,10, 11 and 14 
do approximate the measured isotherms. The fact that the 
predictions are first order estimates with large standard 
errors is immediately apparent. Statistical fitting of data 
which knowingly do not encompass all the determining mechanisms, 
e.g. currents, thermal conductivity of lake, radiation, etc. 
has many shortcomings. In some cases, it would have been 
possible to estimate some of the mechanisms like currents and 
radiation from other surveys or semi-empirical equations. 
However, these estimations were avoided as experience has 
indicated that currents very near to shore and near the surface 
tend to be very local in nature. Particular caution is also 
indicated in the Lakeview area where major shoreline modifications 
like the cooling water intake and wharf structure have taken place, 
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The method is a first order estimate based on 
statistical fitting of limited data taken at one operating 
thermal generating station. It is suggested that the errors 
of estimate should be presented with the estimates to clearly 
identify the limitations of the methods. Transposing the 
predictions to another area has not been tested. 

CONCLUSIONS 

A method for predicting the isotherms of a surface 
cooling water discharge was developed statistically from 
temperature surveys of an existing thermal plume on Lake 
Ontario, However, the errors of estimating the plume were 
large indicating that a more comprehensive measuring program 
is needed to improve the predictions. The measuring program 
should be continuous over many hours and include measurements 
of all the important heat exchange mechanisms. The heat lost 
by mixing with colder lake water was significantly greater than 
the heat lost to the atmosphere (not considering evaporation 
implicitly) . 
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WIND COMPONENT ALONG SHORELINE (Wp) 
RADIUS OF PLUME CENTERLINE (R 
LAKE TEMR (1^) - TI.T'F 



COOLING WATER OUTFALV?*^ 

7.65 m.p.h. 
3,830 FT. 





WIND - 10 m.pji. 



FIGURE 14 
IBSERVED AND PREDICTED (UNIVERSAL EQUATION) ISOTHERMS (JULY 23, 1969) 



TABLE 1 
Statistical Fitting 
Plume Area (A) and Distance from Outfall (X) 



Data 
Set 


Equation 


Std. 

Dev. 

A 




F 


i 

F § 
a=0.01 


Multiple 
Correlation 


Std. 
Error 

1 


Summer ^ 


A» 


= -0.696 + 1.384X* 


4.62 


F 


=117.1 

» 9 5 


6.91 


0.74 


3.11 




A 


= 2.11 + 0.124X^ 


4.62 


F 


= 96.2 

» 9 5 


6.91 


0.71 


3.28 




A 


= 4.38 + 0.00029e^ 


4.62 


F 


= 28.5 

f 9 5 


6.91 


0.48 


4.07 


Winter^ 


A 


= -4.047 + 2.43X 


5.80 


F 


=283.6 

1 6 3 


7.06 


0.90 


2.49 




A 


= 0.91 + 0.239x^ 


5.80 


F 


=283,3 

t 6 3 


7.06 


0.90 


2.49 




A 


= 4.31 + 0.0033e^ 5.80 

1 1 


F 


=114.0 

/ 6 3 


7.06 


0.80 


3.48 



6 2 

Notes: # A (Area) in 10 ft. 



* X (Distance) in 10 ft. 

^ 'Summer' is defined when lake temperature >^45°F 

^ 'Winter' is defined when lake temperature <45°F 



TABLE 2 



Statistical Fitting 



*« 



,** 



Plume Area (A*) vs Temperature Differences AT ^ (Plume-Lake) and AT (Plume-Air) 



1 

Data , r. ^ • 
set i Equation 


Std. 

Dev. 

A 


F 


F § 
a=0.01 


Multiple 
Correlation 


Std. 
Error 


! 

'Summer'' 'a = 105.16 - 10.83 In (Q*AT ) 
+ 0.53 In (QATp^j 

QATpj^ : t = 10.65; t^ ^^=2.384 

Q%A = ^ = ^•°3' ^0.01 =2.384 

■ 


4.676 

- 


1 

F =72.7 

2» 72 


4.91 


0.82 


2.73 


'Winter'' A = 76.77 - 8.45 In (QAT ) 
; +1.25 In (QATp^) 

,QATp^: t = 8.98; t^^^ =2.395 

,QATp^: t = 2.17; t^^^^ =2.395 


5.194 

1 


F =42.6 

2, 57 


5.00 

/ 1 


0.77 


3.35 



NOTES 



* A (Plume Area) x 10 ft^. 

# Q (Discharge from outfall) in cfs 



*# AT (Temperature difference between plume and lake) F 
#* AT (Temperature difference between plume and air) F 
'Summer' is defined when lake temperature >45°F 

2 O 

'Winter' is defined when lake temperature <45 F 



TABLE 3 



Statistical Fitting 

Plume Area (A) and (QAT^^ ) 
(2,500 cfs > Q > 1,400 cfs) 



Data 


fstd. 
Equation Dev. 


F 


1 

F @ 
a= 0.01 


f 

Multiple 
Correlation 


Std. 
Error 


'Summer ' ^ 


A» 


= 14.56 - 1.06 (AT**) 


4.62 


Fi, 95-112. 4 


6.91 


0.74 


3.14 




A 


= 14.14 - 0.00055 
(Q*ATpj^) 


4.62 


Fi,95= 98.4 


6.91 


0.71 


3.26 




A 


= 23.67 - 8.79 In {AT„^ ) 


4.62 


Fi, 95=142. 5 


6.91 


0.77 


2.94 




A 


= 90.53 - 8.86 In (QATp^ 


4.62 


Fi, 95=135. 2 


6.91 


0.77 


2.99 




A 


= 23.53 - 6.28 /Fy^^ 


4.62 


Fi, 95=130.0 


6.91 


0.76 


3.02 




A 


= 23.25 - 0.14 v^Tpj^ 


4.62 


Fi, 95=118. 9 


6.91 


0.75 


3.10 


'Winter " ' 


A 


= 15.91 - 0.89 (ATpj^) 


5.80 


Fi, 63=102. 9 


7.06 


0.79 


3.60 




A 


- 13.97 - 0.00034 (QAT„) 


5.80 


Fi,63= 59.8 


7.06 


0.70 


4.19 




A 


= 23.70 - 7.81 In (AT„,) 


5.80 


Fi ,6 3 = 176.8 


7.06 


0.86 


2.99 




A 


= 74.00 - 6.88 In (QA'^)5.80 


Fi, 63=126. 3 


7.06 


0.82 


3.37 




A 


= 24.76 - 5.81 ^Afpj^ 


5.80 


Fi, 63=155. 4 


7.06 


0.84 


3.14 




A 


= 21.75 - 0.11 /qEh^^ 


5.80 


Fi,63= 94.3 


7.06 


0.77 


3.70 



Notes : 



# 



A (Area) in 10^ ft^ 

Q (Discharge) in cfs 

ATpj^ (Temperature difference between plume and lake) in 



'Summer' is defined when lake temperature 
'Winter' is defined when lake temperature 



>45°F 
<45'*F 



TABLE 4 



Statistical Fitting 
Plume Distance (X) # and Temperature Difference (A": 



PL^ 



Data 
Set 


Equation 


Std. 

Dev. 

X 


F 


F § 

a= 0.01 


Multiple 
Correlation 


Std. 
Error 


'Summer' * 


X = 7.49 - 


0.38 (ATp^) 


2.48 


Fi, 95 = 29. 8 


6.91 


0.49 


2.17 




X = 7.39 - 


0.00020 (Q#*ATpj^) 


2.48 


Fi ,95=28.5 


6.91 


0.48 


2.19 




X = 10.25 - 


2.9 In (ATpj^) 


2.48 


Fi, 95 = 27. 9 


6.91 


0.48 


2.19 




X = 33.08 - 


3.0 In (QATp^) 


2.48 


Fi, 95=29.0 


6.91 


0.48 


2.18 




X = 10.49 - 


2.16 /ATp^ 


2.48 


Fi, 95 = 29. 5 


6.91 


0.49 


2.18 




X = 10.48 - 


0.05 /QATpj^ 


2.48 


Fi, 95 = 29. 5 


6.91 


0.49 


2.18 


•Winter'^ 


X = 7.73 - 


0.32 (ATpj^) 


2.16 


Fi,63= 87.8 


7.06 


0.76 


1.40 




X = 7.00 - 


0.00012 (QATp^) 


2.16 


Fi,63= 51.3 


7.06 


0.67 


1.61 




X = 10.16 - 


2.64 In (ATpj^) 


2.16 


Fi,63= 98.4 


7.06 


0.78 


1,36 




X = 26.88 - 


2.30 In (QATp^) 


2.16 


Fi,63= 73.38 


7.06 


0.73 


1.48 




X = 10.73 - 


2.0 3 /ATpj^ 


2.16 


Fi ,63 = 106.77 


7.06 


0.79 


1.32 




X = 9.61 - 


0.04 /QATp^ 


2.16 


Fi,6 3= 67.78 


7.06 


0.72 


1.51 



Notes 



# 

* 

#* Q 
1 

2 



X (Distance from outfall measured along plume centreline) in 10^ ft 
AT„ (Temperature difference between plume and lake) in °F 

(Discharge from outfall) in cfs. 
'Summer' is defined when lake temperature >.4 5°F 
'Winter' is defined when lake temperature <45°F 



TABLE 5 

Statistical Fitting 

#* 
Discharge/Plume Area (Q*/A#) and Temperature Difference (AT ) 

(For any discharge values) 



Data 
Set 



Equation 



Std. 
Dev. 



F 



F§ Multiple Std. 
a=0.01 Correla- Error 
- tion 



Summer' 



ln(Q/A) = 1.91 + 2.17 InAT 
ln(Q/A) = -12.5 + 2.3 InAT 



PL 



+ 1.9 In Q 



PL 



'Winter' In (Q/A) = 3.11 + 1.42 InAT 

ln(Q/A) = 20.3 + 1.58 InAT 
- 2.3 In Q 



PL 



PL 



1.30 


F 

1 ; 9 5 


^ 


82 


6.91 


0.68 


0.96 


1.30 


F 

*2;95 


° 


46 


4.83 


0.70 


0.93 


1.19 


^1;63 


= 


89 


7.06 


0.76 


0.77 


1.19 


^2;62 


■= 


53 


4.97 


0.79 


0.73 



Notes : 



# 

* 



1 
2 



A 
Q 



(Area) in 10^ ft^ 
(Discharge) in cfs 



PL 



(Temperature difference between plume and lake) in °F 



'Summer' is defined when lake temperature M5°F 
'Winter' is defined when lake temperature <45 °F 



APPENDIX I 
Sample Computations for Thermal Plume Predictions July 23, 1969 

Relevant Data 

Air temperature, T - 75°F 
Lake water temperature, T - 71. 7° F 
Wind Speed, W = 10 mph @ 2 71° 
Outfall Discharge, Q = 1750 cfs 

Definition of Plume Centreline 

Since the shoreline is @ 231° , the angle 9 between 
shoreline and prevailing wind is: 

= 271-231 = 40° 
. . Wind Component along shoreline 

Wp = W Cos = 10 Cos 40 = 7.6 5 mph 
But R = 500 |Wp| ft 
. . R = 3830 ft. 

The arc of the plume centreline is drawn with a radius 
of 3830 ft. on locii of centres (see Figure 4). Since the 
arc does not lie wholly on the water, the plume centreline is 
smoothly joined to outfall, avoiding the shore. 

Definition of Plume Area 

Since the lake temperature is 71. 7° F, the computations 
are done for 'suiraner' regime. The prediction equation is 
accordingly selected from Table 1. 

A = -0.696 + 1.38 4X (1) 

6 2 
where A = area in 10 ft. 

X = distance from outfall, measured along 

-3 
the plume centreline in 10 ft. 

The following computations are for X = 4 (4000 ft.) 

(a) Mean Area 

A = -0.696 + 1.384(4) = 4.84 



(b) Maximum Area (One standard error) 

Standard error of equation (IJ = 3.11 from Table 1. 

. . Maximum area A* = 4.84 + 3.11 = 7.95 

Determination of Temperature Difference ("^Tpj^) 

The prediction equation is selected from Table 3 as 
follows : 

A = 90.53 - 8.86 In (Q^Tpj^) (2) 

where ATp^ = T^ - T^ in °F 

T = plume temperature in °F 

(a) '^T (Mean Area) 



Substituting the value A = 4.84, calculated above 
for X = 4 and Q = 1750 in equation (2) , ^t^^ is computed. 

. . ATpj^ =9.05 

.'. T^ = T, + At„^ = 71.7 + 9.05 = 80.75°F 
P L rL 

(b) At * (Maximum Area) 

Using A* = 7.95 and Q = 1750 in equation (2), 
At„ * is computed. 

XT Lj 

At * = 6.37°F 

•. T * = T + At^ * = 71.7 + 6.37 = 78.07°F 
P L PL 

Determination of Isothermal Contour 

Equation (3) describes a multi-leaved rose in polar 
co-ordinates (r,ct ) 

r = x Cos n" • 

where n governs the number of leaves and therefore the 
lateral spread. (See Figure 1) . 

Area enclosed by the curve may be shown as : 

A = :ix! <^^ 

4n 

For X = 4, substituting the values of A = 4.84 and 
A* = 7.95 in equation (3) values of n and n* are computed 

n = 2.6 

n*= 1.58 
Substituting these values of n and n* in equation (3) , the 

corresponding curves can be plotted. 



(3) 



APPENDIX II 
SAMPLE COMPUTATIONS FOR THERMAL PLUME PREDICTIONS - JULY 23, 1969 

Relevant Data 

Air temperature, T = 75°F 

Lake water temperature, T^ = 71.7°F 

Wind speed, W = 10 mph @ 271° 

Outfall Discharge, Q = 1,750 cfs 

Definition of Plume Centreline 



Since the shoreline is ^ 231° , the angle 9 between shoreline 
and prevailing wind is: 

e = 271-231 = 40° 

.*. Wind component along shoreline 

W = W Cos e = 10 Cos 40 = 7.65 mph 

P 

But R = 500 |Wp| ft 

. * . R = 3,830 ft. 

The arc of the plume centreline is drawn with a radius of 3,8 30 

ft on locii of centres (see Figure 4). Since the arc does not 

lie wholly on the water, the plume centreline is smoothly joined 

to outfall, avoiding the shore. 

Definition of Plume Area 

Since the lake temperature is 71.7°F, the computations are done 
for 'summer' regime. The prediction equation is accordingly 
selected from Table 1. 



A = -0.696 + 1.384X • • • (D 

where A = area in 10^ f t ^ 

X = distance from outfall, measured along the plume 
centreline in 10 ^ ft. 
The following computations are for X = 4 (4,000 ft) 

(a) Mean Area 

A = -0.696 + 1.384 (4) = 4.84 

(b) Maximum Area (One standard error) 

Standard error of equation (1) = 3.11 from Table 1. 
.*. Maximum area A* = 4.84 + 3.11 = 7.95. 

Determination of Temperature Difference (AT p^) 

The prediction equation is selected from Table 5 as follows 

ln(Q/A) = 1.91 + 2.17 In(ATpj^) ... (2) 

where AT_,, = T^, - T^ in °F 
PL P L 

T = plume temperature in °F 

(a) AT^, (Mean Area) 

Substituting the value A = 4.84, calculated above for X = 4 
and Q = 1,750 in equation (2), ATp^^ is computed. 

.* . ATpj^ = 6.26 

.-. T„ = T, + AT„, = 71.7 + 6.26 = 77.96 " 78.0°F 
P Ij PL 



APPENDIX III 
Data Used 



Survey Date Discharge Outfall Lake Air 

Flow Temp. Temp. Temp, 



J? O— o», o_ 
cfs F F F 



Apri 


1 4, 1969 


1970 


56.2 


40.0 


52.0 


May 


22, 1969 


1400 


62.0 


46.0 


50.0 


May 


23, 1969 


1400 


61.0 


46.0 


57.0 


June 


19, 1969 


1830 


58.6 


47.0 


60.0 


July 


15, 1969 


17 30 


76.0 


61.7 


82.0 


July 


16, 1969 


1730 


78.0 


62.2 


83.0 


July 


17, 1969 


1620 


74.0 


61.7 


83.0 


July 


22, 1969 


1750 


84.0 


69.3 


75.0 


July 


23, 1969 


1750 


84.0 


71.7 


75.0 


July 


24, 1969 


1750 


86.0 


71.5 


76.0 


Sept 


, 9, 1969 


1970 


72.0 


59.0 


66.0 


Sept 


. 10, 1969 


1750 


70.0 


58.0 


69.0 


Sept 


11, 1969 


1750 


74.0 


59.0 


63.0 


Sept. 


30, 1969 


1830 


73.0 


58.0 


66.0 


Oct. 


1, 1969 


1830 


73.0 


5 8.0 


60.0 


Oct, 


2, 1969 


1980 


75.0 


59.5 


61.0 


Oct. 


21, 1969 


1740 


62.0 


43.0 


48.0 


Oct. 


22, 1969 


1760 


58.0 


43.0 


31.0 


Oct. 


23, 1969 


1720 


56.0 


42.0 


33.0 


Nov. 


12, 1969 


2130 


64.0 


47.0 


47.5 


Nov. 


13, 1969 


2140 


60.0 


46.5 


42.5 


Nov. 


25, 1969 


2040 


57.0 


41.0 


45.0 



(b ) Tp (Maximum Area) 

Using A* = 7.95 and Q = 1,750 in equation (2), AT_,^ * is 
computed. 

^^PL* " ^'^^ * 5.0°F 



Determination of Isothermal Contour 

Equation (3) describes a multi-leaved rose in polar co-ordinates 

(r,a) 

r = X Cos na , , . (3) 

where n governs the number of leaves and therefore the lateral 

spread (see Figure 1) . 

Area enclosed by the curve may be shown as: 

A = TTXJ_ ... (4) 

4n 
For X = 4, substituting the values of A = 4.84 and A* - 7.95 
in equation (3) values of n and n* are computed n = 2.6 
n*= 1.58 

Substituting these values of n and n* in equation (3) , the 
corresponding curves can be plotted. 



APPENDIX III (cont'd) 
Data Used 



Survey Date 


Discharge 
Flow 


Outfall 
Temp. 


Lake 
Temp. 


Air 
Temp. 






cfs 


°F 


Op 


Op 


Nov. 26, 


1969 


2030 


56.0 


41.0 


40.0 


Nov. 27, 


1969 


2010 


56.0 


41.0 


32.5 


Jan. 28, 


1970 


2520 


56.0 


35.0 


36.0 


Feb. 3, 


1970 


2190 


51.7 


34.0 


8.0 


Feb. 4, 


1970 


2150 


51.4 


33.0 


13.0 


Feb. 20, 


1970 


2205 


61.0 


35.0 


22.0 


May 20, 


1970 


1924 


63.0 


45.0 


52.0 


May 21, 


1970 


1920 


64.0 


46.0 


55.0 


June 10, 


1970 


1860 


66.0 


61.0 


71.0 


June 1 1 , 


1970 


1850 


66.0 


54.0 


71.0 


June 11, 


1970 


1740 


65.0 


54.0 


84.0 


Sept. 1, 


1970 


2270 


65.0 


56.0 


6 5.0 


Sept. 2, 


1970 


2440 


6 8.0 


59.0 


60.0 


Sept. 2, 


1970 


2380 


67.0 


55.0 


65.0 


Oct. 20, 


1970 


2080 


70.0 


54.0 


53.0 


Oct. 20, 


1970 


1980 


72.0 


54.0 


54.0 


Oct. 22, 


1970 


1970 


72.0 


55.0 


55.0 


Oct. 22, 


1970 


1630 


72.0 


55.0 


56.0 
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